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We investigate one-electron properties of one-dimensional self-similar structures called limit quasi- 
periodic lattices. The trace map of such a lattice is nonconservative in contrast to the quasi-periodic 
case, and we can determine the structure of its attractor. It allows us to obtain the three new features 
of the present system: 1) The multi-fractal characters of the energy spectra are universal. 2) The 
supports of the /(a)-spectra extend over the whole unit interval, [0, 1]. 3) There exist marginal 
critical states. 
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Since the discovery of a quasicrystal, deterministic and 
aperiodic systems have been attracted much attention 
p). They are classified into the third group in addi- 
tion to periodic systems and random systems. Although 
they include a wide range of diverse structures, a self- 
similar structure such as the Fibonacci lattice or the Pen- 
rose tiling is especially important because quasicrystals 
belong to materials of this type. In particular, a one- 
dimensional (ID) self-similar structure has the symmetry 
described by a semi-group,and the electronic state on it 
exhibits rich properties yl El The presence of mul- 
tifractal energy spectra and self-similar wave functions is 
peculiar to such systems. Here we will add universality 
to the subject; it applies to limit quasiperiodic lattices 
which form an important class of self-similar structures. 

A ID self-similar structure is produced with a substitu- 
tion rule (SR) starting from two elements, A and B 
It is taken to be an aperiodic lattice composed of two 
types of sites, and called a self-similar lattice (SSL). A 
general form of the SR is written as A — > <ta(A, B), B — > 
ctb(A, B), where <ja{A, B) and erg (A, B) are words of 
the two letters or, equivalently, finite sequences of the 
two symbols. The SR can be specified by a combina- 
tion, (<ta(A, B),ctb{A, B)), of two different words, and 
an infinite number of different SRs are formed by chang- 
ing the combination. A simplest SSL is the Fibonacci 
lattice, whose SR is (B, AB), while another one with 
the SR, (B, BAAB), produces the mixed mean lattice 
(MM lattice) Q- To every SR, (a a, <Jb), we can as- 
sociate a Frobenius matrix, M :— {^ a c b A with a and 

c (resp. b and d) being the numbers of A and B in 
(ta{A, B) (resp. ctb(A, B)); we shall called it the sub- 
stitution matrix of the SR or of the relevant SSL. An 
SR is not uniquely determined by its substitution ma- 
trix because two or more choices are usually possible for 
the orders of the letters, A and B, in a a and erg. If 
there are several locally-isomorphic classes of SSLs with 
a common substitution matrix, we may call them iso- 
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mers. For example, two SSLs with SRs, (B,ABBA) and 
(B,ABAB), are isomers of the MM lattice. When an 
SSL has a global center of the reflection symmetry, it can 
be called a symmetric SSL and the relevant SR a sym- 
metric SR 0. We may call an SR, (a a, <Jb), palindromic 
if both a a and ob are palindromes, while we may call it 
quasi-palindromic if it is written with a palindromic SR, 
(a' A , a' B ), as (a' A B, <j' b B). Palindromic SRs and quasi- 
palindromic ones produce SSLs. If a substitution ma- 
trix is specified, at least one of the relevant isomers is 
a symmetric SSL. All the SRs presented above produce 
symmetric SSLs. 

We shall turn to the structure factor S(Q) an SSL 
0,^3- The set of all the SSLs is divided into two classes 
LT and T; S(Q) for the case IT consists only of Bragg 
peaks, while the one for the case T is singular continu- 
ous. The class II is divided, further, into the three sub- 
classes, i) LIi : Quasiperiodic, ii) Iln: Limit Quasiperi- 
odic, and hi) Hn: Limit Periodic. Representatives of 
the three are the Fibonacci lattice, the MM lattice, and 
the period-doubling lattice, respectively. Let M be the 
substitution matrix of an SSL, C, and r (> 1) be the the 
Frobenius eigenvalue of M. Then, a necessary and suffi- 
cient condition for C to belong to LTi or Iln is that the 
second eigenvalue r' satisfies |r'| < 1. If it is satisfied, 
r is a quadratic irrational called a Pisot number. Note 
that t' is the algebraic conjugate of r and rr' = s with 
s := detM; M is unimodular only when s = ±1. It is 
known that £ belongs to LTi if |s| = 1 but to Iln other- 
wise. For example, r = (1 + \/5)/2 (resp. r = 1 + \/3) 
for the Fibonacci lattice (resp. the MM lattice). The 
set of all the wave numbers at the peaks of S(Q) form 
an additive group, M., which is called the Fourier mod- 
ule. The two subclasses Ii\ and Iln are distinguished by 
whether the number of the generators of M. is two or in- 
finite. Let M be the substitution matrix of an SSL and 
let lj := (t — a)/b — c/(r — d) with r being the Frobenius 
eigenvalue of M. Then, the Frobenius eigenvector of M is 
given by the column vector *(1, oS). The Fourier module 
M of an SSL in Iln is written as M = (1 + uj)~ 1 Z{oj} 
with Z{lu} := Z[w] U t~ 1 Z[lu] U t~ 2 Z[lu] U • • • , where 
Z[lu] := {x+yuj \ x,y S Z}. The Fourier module is a dense 
set on the real axis. Note that M. is common among dif- 
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FIG. 1: The energy spectrum of the MM lattice with A := 
Vb-Va = 1 



ferent isomers because M. is uniquely determined by M. 
Incidentally, Z{uj} = Z[w] and M = (1 +w)~ 1 Z[cj] for 
an SSL in IIi becasue rZ[w] = Z[w] for this case. 

Every set, say II, of SSLs is divided into the symmetric 
subset II s and its complement IP, where II s is composed 
only of symmetric SSLs. Previous researches on the one- 
electron states of SSLs are focused exclusively on the case 
of class IP. Moreprecisely, the subclass ITf is investi- 
gated extensively [5| but Ilfj is investigated by nobody, 
while there are a considerable number of papers on the 
subclass IIjjj but they are still on a rudimentary stage 
(see, for example, 0). Therefore, the authors investi- 
gated in detail the case of the subclass Ilfj. Though ITf x 
appears a simple extension of II| , it turns out as shown 
below that its one-electron states exhibit a remarkably 
different character from the case of the latter. We shall 
call hereafter an SSL in Ilfj a limit quasi-periodic lattice 
(LQL). 

Self-similarity of an SSL allows us to use the trace map 
(TM) for researches of the electronic states on it 0, H|. 
The character of the TM as a non-linear ma p depe nds 
crucially on whether it is conservative or not ^ 
is known that SSLs with conservative TMs form a special 
subset, nj, of nf 0,112. Hereafter, an SSL in Uf is called 
simply a conservative SSL (CSSL). Remember that an 
LQL is never conservative. 

Generic characters of one-electron states of an SSL can 
be investigated using the tight-binding Hamiltonian 7i, 
each site-energy of which takes Va or V B depending on 
the type of the site [5| • The unit of energy can be so cho- 
sen that the transfer integral assumes — 1. It was shown 
generally that the energy spectrum <r = cr(H) is singu- 
lar continuous 0] . More precisely it is a bounded closed 
set having an infinite number of gaps, and, moreover, the 
length of its every connected component vanishes; we can 
see it in Fig. 1 . This means that the integrated state den- 
sity H{E) per one site is a devil's staircase. Every gap 
of a can be specified by the height of the corresponding 
step of the staircase. let G a be the set formed by adding 
and 1 to the set of the heights of all the steps. Then, 
it is a discrete subset of the unit interval [0, 1], and is 
determined by the gap labeling theorem 14] . Since a is 
a multifractal, its singularity, a — a(E) with E G er, is a 
function of E. The distribution of a in a is characterized 
by the /(a)-spectrum The support of f(a) is a closed 
interval, [a m i n ,a max ] C [0, 1]. For every E in a, we can 
define a function by h(e) := H(E + e/2)- H(E-e/2). If 



it satisfies the equation, h(s) ~ fih(e/X), for two param- 
eters A = X(E) and \i = n(E) chosen suitably, E is a cen- 
ter of local self-similarity of a. The singularity a of the 
center is given by a — In /i/ In A 3, 5j. A remarkable fea- 
ture of an SSL is that the set, a\ s , of all such centers is a 
countably infinite set 0. The set, E := {H(E) | E e ai s } 
is an important object which characterizes a\ s . 

In the trace map formalism, a homomorphism is de- 
fined from J- (A, B), the free semi-group generated by 
the two generators A and B, into SL(2, R), a group 
formed of 2 x 2 unimodular matrices 010. Ifi 
and B are the images of A and B in the homomor- 
phism, the image, B), of a sequence W(A, B) is 
called the transfer matrix associated with the sequence. 
We can form an infinite series of pairs, {A n , B n }, of se- 
quences by using the recursion relation {A n+ i, B n+1 } — 
{<JA(A n , B n ),<TB(A n , B n )} supplemented by the initial 
condition, {Ai, Bi} := {A, B}. The corresponding pair, 
{A n , B n }, of transfer matrices satisfy a homologous re- 
cursion relation. Let {A, B} be such a pair at any gener- 
ation and W(A, B) any sequence. Then, TrH 7 ^, B) is 
a polynomial with integral coefficients of the three vari- 
ables, x := TrA, y := TrB, and z := Tr(AB). Thus, 
a map is defined from J- (A, B) into Z[x, y, z]. Let T x , 
T y , and T z be the three polynomials defined as the im- 
ages of the three sequences, <ta{A 7 B), ctb(A, B), and 
(ja(A, B)<tb(A, B), respectively. Then, the trace map 
T which yields a 3D dynamical system over a 3D phase 
space R 3 is defined byr -> Tr := (T x (r), T„(r), T z (r)) 
with r = (x, y, z) 0,|E EI- For example, Tr assumes 
(y, x, yz — x) for the Fibonacci lattice but 

(y, xyz - x 2 - y 2 + 2, x{y 2 - l)z ~ x 2 y - (y 2 - 3)y) (1) 

for the MM lattice. Noncommutativity between the 
transfer matrices in a pair JA, B} can be quantified by 
the polynomial defined by 0] 

I := -Tr[A, B} 2 = x 2 + y 2 + z 2 - xyz - 4, (2) 

which vanishes if they are commutative. It can be shown 
readily that J(ri) = (A/2) 2 . If the subscript n be- 
ing suppressed at the moment is recovered, we may 
write r„ + i = Tr n . An orbit 0(E) := {r n \n e N} 
in the phase space is generated from the initial point 
r x = (V A -E,V b - E, (V a - E)(V B - E) - 2). An im- 
portant property of the TM is that Tr itself does not 
depend on the parameters Va, Vb, and E. 

It can be shown generally that there exists a polyno- 
mial P(r) such that [A n+1 , B n+1 ] = P(r n )[A n , B n ] for 
a palindromic SR but [A n+1 , B n+1 ] = P(r n )[A n , B n ]B n 
for a quasi-palindromic one f!7| . Moreover, |P(r)| = 1 
for the conservative case but deg(P) > 1 otherwise. We 
obtain P = z for the MM lattice, for example. It yields 
the identity /(Tr) = [P(r)] 2 /(r). Therefore, I(r) is an 
invariant of the TM only for the case of a CSSL but I(r) 
is a semi-invariant for the general case because the TM 
never changes its sign [||. 
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Let us assume that A n and B n for 3 neN are not com- 
mutative and, besides, that P(x n ,y n , z n ) — for some 
E 6 (7. Then, A m , B m are commutative for v m > n + 1, 
and we can conclude that the relevant eigenstate will be 
extended as in the case of a periodic system 0] . Differ- 
ent extended states may be obtained from the roots of 
a similar equation at a different generation. Therefore, 
an SSL may have an infinite number of extended states 
unless it is not an CSSL. 

The dynamical system defined by the TM, T, is charac- 
terized by its limit cycles (LCs). If a point in R 3 belongs 
to an LC with period p, it is a fixed point of T p and an 
orbit starting from it is a pure cycle. A necessary and 
sufficient (N&S) condition for E 6 a to belong to <j\ s is 
that 0(E) falls on an LC. If this is satisfied, we can de- 
termine the scaling parameter A by a linear analysis of 
T p around the relevant fixed point, where p is the pe- 
riod of the LC, so that the singularity a = a(E) can be 
evaluated The wave function of the corresponding 
eigenstate is known to be asymptotically self-similar; the 
spatial ratio of the self-similarity is equal to t p @ . 

Though I(x, y, z) is not an invariant of the TM, T, 
of an LQL, the curved surface defined by the equation 
I(x, y, z) = is an invariant surface. Its restriction, S, 
into the cube [—2, 2] 3 e R 3 is a closed surface, which 
is a ball with the tetrahedral point symmetry but has 
four cusps at (2, 2, 2), (2,-2, -2), (-2, 2, -2), and 
(—2, —2, 2) (see Fig.lc in [ill). A coordinate system can 
be introduced into S by x — 2 cos (2ttu), y = 2 cos (2ttv), 
and z — 2 cos [2n(u + v)] with with (u, v) being a vari- 
able on the 2D torus T 2 := R 2 /Z 2 . Actually, S is doubly 
covered by T 2 because ±(u, v) are mapped onto a sin- 
gle point on S. The origin (0, 0) of T 2 corresponds to 
(2, 2, 2) in R 3 . The TM induces a 2D dynamical sys- 
tem (DS) over T 2 : T(u, v) = (u, v)M with M being 
the relevant substitution matrix. With these machiner- 
ies together with a theory of algebraic number theory, we 
can determine all the cycles of the 2D DS, which will be 
presented elsewhere. 

For a periodic system (A = 0), every 3D orbit of T 
is confined to S from the very beginning. The initial 
state ri of the orbit 0(E) with E = — 2cos27rre corre- 
sponds to (ui, Vi) = ±(re, k) on S, where re is the ra- 
tionalized wave number of a plane wave state. It follows 
that (u n , v n ) = ±(L„ k, Ln*^ k) with Ln an d in be- 
ing the numbers of the letters in A n and B n , respectively. 
A N&S condition for the 2D orbit {(«„, v n ) \ n <E N} to 
converge on (0, 0) € S is equivalent to re G M 0,113. On 
the other hand, the 2D orbit falls on a limit cycle with 
period p if and only if {(u n+p =p u n , v n+p =p v n ) \ n e N} 
converges on (0, 0) € S. Using the fact that the ra- 
tios L^p/L^' and L^ p /L^ tend to t p as n goes 
to the infinity, we can show that a N&S condition for 
that condition to be satisfied is that (t p — l)re € M. or 
(t p + 1)k € M. It follows that a N&S condition for the 
2D orbit to fall on a limit cycle is given by re £ Q [w] — M. 
with Q[w] := {u + vlu \ u, v G Q}. Thus, the behav- 
ior of the TM of a periodic system has been completely 



revealed. Since the energy spectrum of the periodic sys- 
tem is absolutely continuous, a(E) takes 1 15]. Two 
states with wave numbers ire degenerate in energy and 
the corresponding two orbits of the TM are identical, so 
that we can assume < re < 1/2. The number re is noth- 
ing but the normalized value per one site of the number 
of the nodes in the relevant sinusoidal wave. Since H(E) 
is related to re = k(E) by the equation H — 2k(E), the 
asymptotic behavior of the TM at the energy E is de- 
termined by the value of H(E), which belongs to the 
interval [0, 1]. 

We will proceed to the case of a non-periodic system 
(A / 0). Since an LQL has the reflection symmetry, ev- 
ery energy level has its own parity with respect to the 
center of symmetry. If A is changed adiabatically from 
zero to the present value, every energy level will change 
continuously because the degeneracies being present at 
A = does not matter if the parity is specified. It fol- 
lows that re remains to be the normalized number of the 
nodes. Then, the asymptotic behavior of the TM will be 
determined by the value H — H(E), and the rule is the 
same as what is described above for the periodic system. 
Thus, we can conclude that E = Q[w] n \0, 1], which is 
a similar result to the one obtained in Q for the case 
of the Fibonacci lattice. A remarkable property of the 
TM of an LQL is its universality. It is based on the two 
points: 1) The TM does not explicitly depend on the pa- 
rameters A and E. 2) It is essentially of 3D character, 
and not confined to a 2D manifold as in the case of a 
CSSL. Thus, the attractor is independent of the value of 
A, and identical to that of the periodic system. It follows 
that the /(a)-spectrum of a is universal. Moreover, we 
can conclude that a = 1 for every center of local self- 
similarity of g 15]. For example, a 3D orbit for the MM 
lattice converges to a 1-cycle if re = 1/3 , and the scaling 
parameter A determined by a linear analysis of T around 
the relevant fixed point, (—1, —1, —1), is found exactly 
equal to r, so that a = 1. As a consequence, we obtain 

C^max 1 • 

From a brief perturbative consideration, we may ex- 
pect that G a = M. fl JO, 1], which is guaranteed by the 
gap-labeling theorem [Tj]. The gap labeled by H e Go- 
is bounded by two energy levels located at band edges. 
However, the corresponding fixed point (2, 2, 2) of the 
TM is a singularity of S, and we cannot conclude that 
a = 1/2 for these two energy levels. Surprisingly, these 
states are actually marginal critical states, which we can 
prove by a similar technique to the one used in |l9l ]. It 
follows that Q m i n = 0. This and the result a max = 1 
are consistent with that the /(a)-spectrum is universal 
because there is no reason why the value of a m i n or a max 
must take an odd value. The broadness of the support 
of the /(a)-spectrum means that the energy spectrum 
of an LQL is more inhomogeneous than that of a CSSL, 
which one can feel if he compares Fig.l with the energy 
spectrum of, say, the Fibonacci lattice We show in 
Fig. 2 the /(a)-spectrum of the energy spectrum in Fig.l. 

If there are several isomers, their TMs are different. 
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FIG. 2: The /(a)-spectra of the energy spectra of periodic 
approximants for the MM lattice. The data from the 4-th to 
the 9-th generations are used. The data in the region a < .2 
are cut because their accuracies are insufficient on account of 
the presence of marginal critical states. The inset is a log-log 
plot of 8 n := Qmax — 1 versus the generation number n. The 
asymptotic linearity of the plot concludes that S n tends to 
zero in an inverse power law for n as n goes to the infinity. 
The /(a)-spectra of the Fibonacci lattice is superposed for a 
comparison. 



Their /(a)-spectra belong, however, to a common uni- 
versality class because each universality class of the /(a)- 



spectrum of the LQL can be specified by the substitution 
matrix. 

The TM derived from a CSSL is a 2D map on the 
curved surface I(x, y, z) = (A/2) 2 . Since this surface 
depends on the value of A, the /(a)-spectrum as well 
as a m ; n and a max are not universal. This is consistent 
with the known fact that a m i n > and a max < 1 for 
every CSSL (see Fig.2) 0. The SSLs in the set ITf are 
actually minorities in LT|, and one-electron properties of 
SSLs in llj — llj has never been investigeted yet. The 
present theory is basically applicable to this case; the 
only modification needed is to replace the Fourier module 
by M. = (1 + a;) _1 Z[w]. Thus, among isomers belonging 
to ITf, only the CSSL member is nonuniversal. Inciden- 
tally, the SR, (ABA, ABABA), produces a CSSL but its 
isomer, (ABA, BAAAB), does not. 

Although no substance which takes a structure repre- 
sented by an LQL as its thermodynamically stable phase 
is known, such a structure is realizable as an artificial 
super-lattice. Since a transfer matrix can be used also 
to a ID Schrodinger equation, the results of this letter 
apply to the electronic state in such an artificial super- 
lattice as wellj^. Furthermore, it is basically applicable 
to a propagation in a stratified substance of an ultrasonic 
wave, an electromagnetic wave (light), and a spin wave 
as well. If, for example, a lower marginal critical state 
is used as a channel of light, the speed of light will fall 
dramatically. 

In conclusion, one-electron properties of limit quasi- 
periodic lattices exhibits a universality in contrast to the 
case of quasi-periodic ones, and a universality class is 
specified solely by the substitution matrix. 

One of the authors (K. N.) is grateful to Y. Morita. 
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